QCD sum rules for the pseudoscalar decay constants 
— To constrain the strange quark mass 



Hungchong Kim * 

Institute of Physics and Applied Physics, Yonsei University, Seoul 120-749, Korea 



Abstract 



We study the higher order corrections of quark masses to the Gell- 
Mann— Oakes— Renner (GOR) relation by constructing QCD sum rules ex- 
clusively for pseudoscalar mesons from the axial- vector correlation function, 
i J d^x e:'P-''{0\T[Af,{x)A^{0)]\0). To project out the pseudoscalar meson con- 
tributions, we apply —p^p'^/p^ to this correlation function and construct sum 
rules for the decay constants of pseudoscalar mesons, /tt, fk and /^g. The OPE 
is proportional to quark masses due to PCAC. To leading order in quark mass, 
each sum rule reproduces the corresponding GOR relation. For kaon and r/g, 
the deviation from the GOR relation due to higher orders in quark mass is 
found to be substantial. But the deviation gives better agreements with the 
phenomenology. Our sum rule provides a sensitive relation between //^ and 
rris, which stringently constrains the value for m^. To reproduce the experi- 
mental value for fx, ms is found to be 186 MeV at 1 GeV scale. The /^g sum 
rule also supports this finding. 
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I. INTRODUCTION 



According to PCAC the decay constant of a pion measures the couphng strength 
of the corresponding axial-vector current to the pion. One useful constraint which can be 
derived by combining PCAC with the soft-pion theorem is the Gell-Mann— Oakes— Renner 
(GOR) relation 

- Am,{qq) = mlfl , = 131 MeV) . (1) 

As the pion decay constant /^r is well-known [Q, this GOR relation can be used to con- 
strain either the quark condensate or the quark mass. In particular, the small quark mass 
nig = [rriu + 'md)/2 ~ 7 MeV at 1 GeV scale restricts the quark condensate to be 
{qq) ~ —(225 MeV)'^, the value being widely used in practical QCD sum rule calculations. 
Corrections to the GOR relation are an order of 0{m?^ or higher, which should be small. 

When a strange quark is involved, the corresponding GOR relation may not be a reliable 
constraint. For kaon, one may start from the kaon PCAC relation 

(0|s7^75M|ir+(p)) = ip^fK . (2) 

By taking the divergence and using the soft-kaon theorem, one arrives at the kaon GOR 
relation 

- {mq + ms)[{qq) + (ss)] = m%fl . (3) 

Due to the large strange quark mass m^, one may expect nonnegligible corrections from 
(9(m^) or higher order, and this kaon GOR relation, though exact to leading order in quark 
mass, may not be a good constraining equation for and fx- Indeed, using the conventional 
QCD parameters [for example, see Ref. 0.] 

m, = 150 MeV ; (ss) = 0.8(gg) , (4) 

Eq.(0) yields /i^: = 116 MeV, which is much smaller than its experimental value /^^* ~ 160 
MeV 0. Much smaller rus from the baryon sum rules leads to much larger discrepancy. 
To reproduce /^^* from this kaon GOR relation, the strange quark mass (assuming that 
(ss) is fixed as above) needs to be around 300 MeV, which is evidently too large for the 
current mass of strange quark. Consideration for the r^g case leads to a similar (even larger) 
discrepancy from the r^g GOR relation. 

Then one may naturally attempt to improve the GOR relations especially when strange 
quark is involved. We look for the improvement that provides a reliable relation between fx 
and rris. This improvement of course should not affect the well-established GOR relation in 
u, d quark sector, i.e. Eq.(0). We then look for the tUs in this improvement that reproduces 
/^^*. For this prediction to be reliable, this must be tested in the other case such as r]^. 

As the discrepancy mentioned above comes from the large strange quark mass, it may 
be natural to calculate higher quark mass corrections to the GOR relations. For this pur- 
pose, we use QCD sum rules [^,|| for the axial-vector correlation function. Specifically, we 
wish to construct a sum rule whose phenomenological side contains only the pseudoscalar 
meson contributions while the QCD side is proportional to the quark masses. To leading 
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order in quark masses, we want it to reproduce the GOR relation precisely so that we can 
systematically study the corrections from higher orders in quark mass to the GOR relation. 
Similar QCD sum rule calculations of the kaon decay constant fx did not quite look into 
this higher order corrections and our approach in this work is different from them. 

To illustrate our approach briefly, we start by noting that the correlation function of the 
axial-vector current for example = ^7/^756? contains the two invariant functions, IIi and 
112, defined by 



Il^. = tJ d'xe'P-^0\T[A^{x) A,(0)]|0) , 

= -g^i^Ui + Pf,pull2 ■ (5) 

This is in contrast to the vector correlation function which contains (under the isospin 
symmetry) one invariant function due to the current conservation. PCAC tells us that pion 
and its higher resonances contribute only to 112 while the axial-vector meson contributions 
are contained in both IIi and 112. PCAC further imposes that p^H^^y picks up contributions 
only from the pionic resonances. Therefore, to investigate the properties of a pion without 
having contamination from the Oi meson, it may be useful to consider the correlation function 
defined by 

n(p2)^_^n^^ = ni-p2n2. (6) 

This projected correlation function in the QCD side must be proportional to quark masses. 
Therefore, the sum rule from Il{p^) may provide a useful constraint which directly relates 
quark masses to the pion decay constant. By generalizing this sum rule to the strange 
quark sector, we investigate the quark mass corrections of the order (9(m^) or higher to the 
corresponding GOR relation. This may give stringent constraints for m^. 

This paper is organized as follows. In Section ||, we present the OPE calculation for the 
projected correlation function of the general axial-vector current ^i7^75(3'2- Using this OPE 



up to dimension 6, we construct QCD sum rules for in Section |ITI| and estimate how large 
the higher order corrections in quark mass to the GOR relation. In Section IV, we apply 
this framework to the kaon sum rule. We use this kaon sum rule to determine the strange 
quark mass by using /^^* as an input. In Section 0, we test this in the rjg sum rule 
and see if it reproduces /^g consistent with the phenomenology. We summarize in Section |V| 



II. THE PROJECTED SUM RULE FOR THE GENERAL AXIAL- VECTOR 

CURRENT 

In this section, we calculate the operator product expansion (OPE) for the projected 
correlation function of the general axial-vector current ^i7/x75Q'2 

- ^"^■^ I I d^x e^-^ (0|T[gi(x)7^75g2(x) g2(0)7.75gi(0)]|0) 



p2 



^ / rf^xe^^'-Tr [in^iS^ix, 0)|^75^5i(0, x)] , (7) 
p J 
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where qj{j = 1,2) = u,d ot s, and iSj{x,0) is the quark propagator of the flavor qj in 
the fixed-point gauge. Once the OPE of this correlation function is calculated, it will be 
straightforward to obtain the OPE for pion, kaon and rj^. The advantage of using this 
projected correlation function is that the axial- vector meson contributions are canceled and 
the continuum threshold starts at a larger scale so that the correlation function is well 
saturated by the low-lying resonance. 

A few technical remarks are in order before we proceed to the detailed OPE calculation. 
As n^(p2) is a scalar function, the OPE involves only even dimensional operators. No odd 
dimensional operators can contribute to the sum rule. Also, as the axial-vector current is 
not conserved by the finite quark mass, all the OPE must be proportional to the quark 
masses. Thus we can focus on the operators that contain quark masses. This sum rule in 
the end will give sensitive constraints relating quark mass to the corresponding meson decay 
constant. 

The OPE of the correlator Eq.(|^ can be calculated by the standard techniques in the 



fixed-point gauge [10|. To include the finite quark mass effectively, we calculate the correlator 



in the momentum space, 

n^(p^) = ^ / ^Tr [^,zS,{k)^,zS,{k - p)] . (8) 

In the coordinate space, one has to expand the propagator in rrij and hence it is difficult to 
keep the quark mass terms to all orders at each dimension. The perturbative part can be 
calculated straightforwardly using the free quark propagator 



^S;-{k)=^p^^ (, = 1,2). (9) 



Putting this into Eq. and using the standard techniques of the dimensional regularization 
and the Feynman parameterization, we readily compute the perturbative part of Eq. (H), 



xln[— m(1 — u)p'^ + ml — u{m\ — m\)\ . (10) 

Note, this is an order of 0{m'j) or higher. For the vector correlation function, the mim2 
term will have the negative sign so that when the two quark flavors are equal (gi = ^2); the 
perturbative part vanishes as can be expected from the current conservation. 

For the nonperturbative part, various even dimensional operators can contribute to the 
sum rule. In this work, we will calculate the OPE up to dimension 6. First we consider 
the case when one quark propagator in Eq. (|) is disconnected to form quark condensate [| 
iS^^{k) J d'^x e*'^'^('?a (3;)^f (0)) while the other is remained to be the free propagator. 
By performing the Tayler expansion around = and the Fourier transformation to the 
momentum space, we obtain the expansion for the quark condensate [^,0, 



In our notation, q, /3, 7 • •• refer Dirac indices and a,b,c ■ ■■ color indices. 
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(11) 



One technical remark [IT^jIH] is that the Wilson coefficient of the condensate {qD^^ ■ ■ ■ D^^q) 
is related to the Wilson coefficient of {qq) via 



and similarly 



C, 



■A 



d 



dk 



d 



dk 



^l■n _ 



Cqq{k) 



(12) 



d 



dk 



d 



dk„ 



(13) 



Using this technique, we calculate the OPE involving the quark condensate up to dimension 
6. Recalling that only even dimensional condensates contribute to the sum rule, the 5°"^ part 
of Eq. ([Tl|) gives nonzero contributions when the quark-mass part of Eq. (^) is taken for the 
remaining propagator. Similarly, the 7"'' part of Eq. (pJ| ) contributes when the ^ part of the 
free propagator is taken. Also, since the OPE must be proportional to quark mass we need 
to extract only the quark mass dependent part from various quark condensates (j = 1,2), 



(qj-fxDe^qj) -gxa^iqjqj) 



{qjD^Dpqj) -> -gafs'^iqjqj) , 



24 



{qm) ■ 



(14) 



In doing so, we obtain the OPE containing the quark condensate up to dimension 6, 



{qiqi) , {q2q2) 



— 777,2 



1. 2 
2^ ^ 



TTln 



— ml 

mi(gigi) 7712(^292) 



(15) 



The corresponding OPE for the vector correlation function is obtained by replacing 777,2 — ^ 
—7772 and {q2q2) —>■ ^(^2Q'2)- The ffist term is the one that leads to the GOR relation at the 
order 0{mj). In our calculation we kept the quark mass in the denominator of the free quark 
propagator. Thus, when we analytically continue to the time-like region, the imaginary part 
picks up the pole at p"^ = m'j. This aspect is slightly different from the usual application 
of QCD sum rules where the pole of the nonperturbative OPE is located at = 0. The 
modification of a sum rule result due to this is marginal in usual practices. However, in our 
sum rule with n(p^), we are looking at small strength both in the phenomenological and 
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QCD side. In this case, this small distinction due to the quark mass in the quark propagator 
could be important especially when strange quark is involved. 

The contribution from the dimension-4 gluon condensate (^^Q^'^ should be vanished 
in our sum rule as it is not proportional to quark mass. The other possible dimension 6 
operators would be mirrij (^^Q"^^ = 1,2) and mj{qjgscr ■ Qqj). These however do not 
contribute to our sum rule. The disappearance of the mj{qjgsCr ■ Qqj) contribution can 
be easily understood by recalling that the GOR relation is exact at the order 0{mj). If 
'^ji^jOsO' ■ Qqj) does contribute, then the GOR relation is no longer valid even at the order 
0{mj). Indeed, one can show by a direct calculation that mj{qjgsa -Qqj) does not contribute 
to our sum rule 0. 

Verifying the disappearance of the mirrij (^^Q'^^ contribution (more precisely its imagi- 
nary part) is technically more involved. For a consistent calculation, one needs to consider 
the quark propagators with one and two gluons attached. Using the quark propagators 
interacting with gluons (see Ref. p] for their explicit expressions), we readily compute this 
contribution and obtain 



TX 



du 



'1 — uS^mi + u^m2 



mi + m2 

Z2 



where L = u{l — u)p^ — m\ + u{m\ — m^) . 



(16) 



Note, the numerator is already dimension 6. Thus, up to dimension 6, we can safely drop the 
quark mass dependence in L. Now the imaginary part of the first term involving (1 — u)'^mi 
is proportional to Jq du 6'{up'^) which of course vanishes. By the same reasoning, the second 
term involving M^m2 also does not pick up the imaginary part. Therefore, the dimension 6 
contribution involving the gluon condensate does not contribute to our sum rule. 

Before closing this section, here we give the imaginary part of the OPE [Eqs. ([TO|) and 
(p!5|)] up to dimension 6 for the correlation function of the general axial-vector current A'^ = 
917^7592, 



-lmn(s) 

TT 



2(mf + ml) [mf — m^)'^ 



+ (mi + 



(mi + m2) 



2 (mf-ml)^ 



+ 2K 



m^ 



hqi)6{s - ml) + {q2q2)S{s - ml) 
mi{qiqi)-^5{s - mj) - m2{q2q2) j^S{s - m^) 



- (mi + m2y 



(17) 



Note, this is symmetric under 1^2. From this, it will be straightforward to obtain the 
OPE for the pion, kaon and t] sum rules to be discussed below. 

Now, the advantage of working with the projected correlation function becomes clear. 
This OPE contains a few sources of uncertainties, quark condensates and quark masses. 



^To check this easily, it is useful to use the quark propagator in the coordinate space as given in 
Eq.(8) and Eq.(13) of the first reference in Ref. [13|. The two contributions to this dimension 6 
condensate indeed cancel each other in our sum rule. 
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Other sources of uncertainty such as quark-gluon mixed condensate, gluon condensate, or 
four-quark operator do not appear in our sum rule, which will certainly reduce the error in 
our prediction. Furthermore, {qq) {q = u, d) can be determined rather accurately from the 
pion GOR relation Eq. (|I]) once is fixed. As we will see in the next section, the higher 
order corrections in quark mass are small in this m, d sector. When strange quark is involved, 
we have only the two additional QCD parameters, and {ss). Thus, as long as it is stable 
under the modeling of higher resonances or Borel mass variation, this sum rule can give a 
stringent constraint for m^. 



III. THE PION SUM RULE 

Having calculated the OPE for the general axial-vector correlator, we can easily obtain 
the pion sum rule. To do this, we set qi = u, q2 = d in Eq. ([I7|) . As the normal GOR relation 
Eq. (|l]) is well established, our improvement including higher orders in quark masses is not 
expected to be large for this pion case. In fact, the small quark masses can guarantee this 
expectation but we will check this for completeness. 

Using the isospin symmetry, m„ = rrid = rriq and (uu) = (dd) = (qq), we obtain the 
spectral density (the imaginary part) of the OPE up to dimension 6 , 



1 3t72^ / 4:771^ 

-lmnr(s) = e{s- Ami) + 4m,(gg)5(s - mj) . (18) 

The phenomenological spectral density has pion and higher resonances of pion, namely, 

-Imnf "'^(s) = -f^ml6{s - ml) + {n' contribution) + ■ ■ ■ . (19) 

71 

Note that there is no ai meson contribution in our projected sum rule. Matching the two 
expressions under the Borel transformation and invoking QCD duality for higher resonances, 
we obtain the pion sum rule 



r^o / 9 r 

/ dse-'^^' -Im nr(s) - nf '^"(s 



^ f^mle-</''' = ^ J'" ds e-^/^^^l - ^ 6(3- Am]) - Amq{qq)e-<"^' . (20) 

The first term (the perturbative part) is an order 0{ml) or higher. To the order 0{mq), we 
recover precisely the GOR relation Eq.(|l]) as expected. Therefore, the perturbative part is 
the correction to the GOR relation. 

We use this sum rule to plot the Borel curve for /^r setting m^r = 138 MeV, mg = 7.2 
MeV The GOR relation gives (assuming = 131 MeV) {qq) = -(225 MeV)^. The 

continuum threshold is taken to be 6*0 = 1.7 corresponding to the vr' mass. Figure shows 
the Borel curve for which is quite stable with respect to the variation of the Borel mass 
M^. Furthermore, the resulting is not sensitive to So- Also shown by the straight line 
is the /,r obtained from the GOR relation. The deviation from GOR relation is about 3 %. 
Thus our sum rule including the higher mass corrections does not change the established 
pion GOR relation. 
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IV. THE KAON SUM RULE TO CONSTRAIN THE STRANGE QUARK MASS 



We now construct a sum rule for kaon. The small deviation from the GOR relation 
observed in the pion sum rule can not be guaranteed in this case because of the large 
strange quark mass. In this sum rule, we set qi = q (here q = u or d) and q2 = s in Eq. (|T7|). 
The similar procedure as before leads to the kaon sum rule 



So 



X 



[rrig + rris) - 



2 K - "^')' 



s 



s — [rriq + ms 



mg{qq)j^e 



(21) 



Again, to leading order in iriq (or m^), we recover the kaon GOR relation Eq. (0). The 
first and third terms are therefore the higher order corrections to the kaon GOR relation. 
To estimate how large the corrections are, we plot for fx with respect to in figure |^. 
In this plot, we use mx = 494 MeV and the continuum threshold is set to So = 2 GeV^ 
corresponding to Kq (1430) Q. For the strange quark condensate, we set (ss) = 0.8{qq) [0. 
As for the strange quark mass, a wide range of its value can be found in literatures [|15|-[20|]. It 
is within the range rris ~ 100 — 200 MeV. To cover this range of rris, we plot the Borel curves 
for rris = 100 (solid), 150 (dashed), 200 (dot-dashed) MeV. The corresponding straight 
lines are the ones from the kaon GOR relation Eq. (^. We see that all three Borel curves 
become quite flat for > 1 GeV^, but the extracted fx is quite different from the fx 
given by the kaon GOR relation, about 20 -25 % level. The corrections shift the Borel 
curves upward and it is clear that the fx extracted from them is closer to its experimental 
value fx^' ~ 160 MeV. Thus our improvement helps to achieve a better agreement with the 
experiment. To reproduce /^^* only from the kaon GOR relation, rris should be abnormally 
large like rris ~ 300 MeV, which is obviously not acceptable. 

Most interesting feature is that the stable Borel curve is shifted noticeably as we change 
rris. It means that, once fx is fixed, this sum rule can give a stringent constraint for rUs or 
vice versa. We therefore look for which reproduces the experimental kaon decay constant 
jexpt ^ -j^gg gjj-^pg g^j^j ^qq^ gj.^ ^gjj fixed by the GOR relation for pion, only uncertainty 
in our sum rule comes from (ss). According to the baryon sum rules ^,14], the strange 
quark condensate at 1 GeV is estimated as 

{ss) = (0.8±0.1)(gg) . 



(22) 



Using this input, we obtain the strange quark mass that reproduces the experimental fx at 
M2 = 1 GeV2, 



^However, the sensitivity to 5*0 is very weak. For Sq = 2.5 GeV^, the result is changed only by 2 
% . 
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= 186.5 ± 8.5 MeV . (23) 

The error comes from the uncertainty in (ss). This is the value at 1 GeV scale because 
other parameters used in extracting this is defined at 1 GeV. The Borel stability in extracting 
this is fairly good. Even if we take into account the additional error of an order 5 MeV from 
the continuum threshold 5*0, the extracted rris from this sum rule is sufficiently accurate. 
Higher dimensional operators starts at dimension 8, which are expected to be small. Our rUg 
is somewhat larger than the conventional value of = 150 MeV or the value from hadronic 
T decay = 159 MeV [|T^]. But it is smaller than the one from Ref. ||TH[ and is somewhat 
consistent with the value from the current algebra ratio = 175 ± 25 MeV or with the 



one in Ref. 1119 



Possible modification to our results can be expected from as correction to the perturba- 
tive part. Note however that the dominant piece to our sum rule is the dimension 4 OPE 
containing quark condensates. Contribution from the perturbative part, which is an order 
(9(m^), is about 20 % level in our kaon sum rule. The correction is expected to be 
proportional to fts/vr times the perturbative part [For example, see the p meson sum rule 
in Ref. 0.]. At M = 1 GeV, is around 0.3 0, which gives the correction of an order 
0.1 to the perturbative part. Thus, unless it involves a large numerical coefficient, the 
correction to our results is expected to be small. Nevertheless, in future, it will be interesting 
to calculate the correction explicitly and make sure this expectation. 



V. THE r/8 SUM RULE 

Having determined rris from the kaon sum rule, we now test whether it gives a consistent 
result for the rjg case by constructing a sum rule for the current 

= ^ (m7/.75M + dl^,15d - 257^75^) . (24) 

We recall that the 773 decay constant is defined by the matrix element 

mllmip)) = ^UP, ■ (25) 

Refs. [0,0 give /r,g ~ 157 — 165 MeV. Using rris from the kaon sum rule, we will see if the 
frjg from our sum rule is consistent within this range. 

Taking the divergence of the axial-vector current and the subsequent use of the soft-r/g 
theorem lead to the following r/g GOR relation 

fvs'^m = -^'^q{^(l) - ^'rrisiss) . (26) 

Here again, we have used the isospin symmetry = = rriq and (uu) = (dd) = (qq). 
The OPE for the correlation function 

H(p2) = J d^xe'P-'^{0\TAl{x)Al{0)\0) (27) 

can be easily obtained from the general formula Eq. (|1^. Assuming that the continuum 
contribution starts at the rj' resonance (i.e. 5*0 = 0.92 GeV^), we easily obtain the rjs sum 
rule, 
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/Vr J irri^ V So 

+ 4 e-/-%/l-M _ lm,i^ss)e-'slM^ . (28) 

vr^ Jiml V s 3 

Again, to leading order in quark mass, this reproduces the r^g GOR relation Eq. (^). The 
perturbative terms that are the second order in quark mass constitute the higher order 
corrections. 

The Borel curve for rj^ is shown by the thick dashed line in Figure |^. In this plot, we use 
ms = 186 MeV obtained from the kaon sum rule and assume that m^g = = 547 MeV. 
To show the sensitivity to the continuum threshold, the curve for So = 1.4 GeV^ is also 
shown by the thin dashed line. Around = 1 GeV^, the two curves differ by 5%. The 
Borel stability is not as good as the kaon case but /^g becomes only 5 % smaller at = 2 
GeV^ than the one at = 1 GeV^. As can be seen, the /^g from this sum rule becomes 
30 % larger than what the r^g GOR relation gives. Therefore, /^g picks up a large correction 
from higher orders in quark mass to the 7]^ GOR relation. From this sum rule, we extract 
/^g = 160 MeV at = 1 GeV^. This is very close to its value in literatures 157 — 165 
MeV PIJ2TII even if we take into account the uncertainties due to the continuum threshold 
or the slight dependence on the Borel mass. Therefore, rris = 186 MeV extracted from the 
kaon sum rule consistently reproduces fx and /^g. 

VI. SUMMARY 

In this work, we have constructed QCD sum rules exclusively for decay constants of the 
pseudoscalar mesons, vr, kaon and r^g. The phenomenological side of our sum rule contains 
only the pseudoscalar mesons while the QCD side is proportional to quark masses. Our sum 
rules give a very stable constraint relating the meson decay constants and the quark masses. 
To leading order in quark masses, our sum rules produce the GOR relations, which allows 
us to study the higher order corrections in quark masses. From the kaon sum rule, we have 
obtained a sensitive constraint relation for J'k and m^, which restricts rris ~ 186 MeV in 
order to reproduce the experimental Jk- Such a rris when applied to the rjg sum rule gives 
frjg consistent with its value in literatures. 
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FIGURES 



FIG. 1. The Borel curve for the pion decay constant. The sohd (straight) hne is the f-^ from 
the GOR relation and the dashed (curved) hne is the /t^ from the Borel sum rule of Eq. (pO|). 
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FIG. 2. The kaon decay constants for given are plotted with respect to the Borel mass. 
For each rrig, we plot the curve from the kaon GOR relation and the one from the kaon sum rule. 
The solid lines are for = 100 MeV, the dashed lines are for = 150 MeV and the dot-dashed 
lines are for m2 = 200 MeV. The parallel straight lines are obtained from the kaon GOR relation 
for given while the corresponding curved lines are from the kaon Borel sum rule Eq. ([2l]). The 
Borel curve (or kaon GOR relation) is shifted noticeably as we change m^. Also, higher order 
corrections of quark mass give 20 ~ 25 % shift from the corresponding GOR value. The fx from 
the Borel curves is closer to its experimental value of 160 MeV. 
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FIG. 3. The 773 decay constant with respect to to the Borel mass. In this plot, we use the best 
fitting nis from the kaon sum rule. The solid straight line is obtained from the rj^ GOR relation 
and the thick dashed curve is from the rjg sum rule (with 5o = 0.92 GeV^) Eq. (]28|) . There is 30 % 
change from the % GOR relation. The thin dashed line is obtained with 5*0 = 1.4 GeV^ to show 
the sensitivity to Sq. 
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